Noise-induced macroscopic bifurcations in populations of globally coupled maps 
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Populations of globally coupled identical maps subject to additive, independent noise are studied 
in the regimes of strong coupling. Contrary to each noisy population element, the mean field 
dynamics undergoes qualitative changes when the noise strength is varied. In the limit of infinite 
population size, these macroscopic bifurcations can be accounted for by a deterministic system, 
where the mean-field, having the same dynamics of each uncoupled element, is coupled with other 
order parameters. Different approximation schemes are proposed for polynomial and exponential 
functions and their validity discussed for logistic and excitable maps. 
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Introduction. A fundamental question that arises when 
physical and biological populations are described by 
means of mathematical models is how the collective dy- 
namics is qualitatively affected by random fluctuations 
acting at the microscopic level. While for systems at 
equilibrium a well established thermodynamic theory ex- 
ists, for out-of-equilibrium phenomena, such as those typ- 
ically described by nonlinear dynamical systems, an over- 
all theoretical view is still lacking. In this Letter we will 
address one aspect of this issue, namely the effect of mi- 
croscopic noise on the mean-field dynamics of large pop- 
ulations of globally and strongly coupled identical maps. 
Systems of globally coupled units are meant to model, 
among others, Josephson junction arrays nl , yeast cells 
in a continuous-flow, stirred tank reactor |2| and neural 
cells , and provide a mean field description of spatially 
extended systems with long enough correlation length. 
The effect of noise on one dynamical system has been 
extensively studied, pointing out phenomena such as 
noise-induced bifurcations, stochastic and coherence res- 
onance (for a review, see Ref. Q and references therein). 
Synchronization phenomena induced by (common or in- 
dependent) noise have been studied for a system of 
two coupled dynamical systems 0, IS 0> IS El E3- I n 
the context of populations, the effect of microscopic 
noise on the collective dynamics has been recently in- 
quired for continuous-time systems, namely phase mod- 
els 0, 0, 0, 113 . integ rate-and-fire neurons flil Hif . 
excitable svstems|rR Il8| . chaotic systems 0, \'2{\ . The 
onset of collective oscillations and the dependence of 
their frequency from the noise intensity has been mainly 
addressed. Analogous synchronization phenomena have 
been detected for stochastic oscillators 0, and in 
spatially extended systems US Ej| . 
Concerning maps, the relationship between single ele- 
ment and mean field fluctuations close to bifurcations [24| 
and the anomalous scaling of the po pul ation moments 
close to the onset of synchronization [2(]j have been in- 
quired. 

However, the qualitative changes of the macroscopic dy- 



namics under the influence of noise have, at our best 
knowledge, never been systematically addressed. In par- 
ticular, the strong coupling regimes could be considered 
of little interest since one expects the addition of a weak 
noise term not to greatly alter the mean field behavior. 
Instead, in the first part of this Letter we give numerical 
evidence of the fact that noise can qualitatively change 
the macroscopic dynamics of the system: even if the sin- 
gle map has a noisy temporal series, the mean field dis- 
plays a low dimensional behaviour which can be different 
with respect to the uncoupled map dynamics. The phe- 
nomenon of noise-induced macroscopic bifurcations will 
be illustrated for two kinds of maps: logistic maps in the 
chaotic regime and "excitable maps" . Contrary to noise- 
induced synchronization, in this case it is the coupling 
that induces a coherent behavior for low noise intensity. 
The deterministic trajectory of each dynamical system 
is blurred out, so that the bifurcations detectable at a 
macroscopic level are a purely collective effect and cannot 
be inferred by looking at one or few population elements. 

In the second part of the Letter, we explain these phe- 
nomena by means of an order parameter expansion, valid 
for sufficiently large coupling strength and large popula- 
tion size. This provides an approximate description of 
the mean field dynamics in terms of few effective macro- 
scopic variables, whose deterministic equations of mo- 
tion account for the macroscopic dynamics of the popu- 
lation and for the bifurcations among different collective 
regimes. The validity of such approach will be demon- 
strated on the two aforementioned populations of maps. 
In particular, for the logistic maps the bifurcation dia- 
gram of the mean field will be rescaled to that of a single 
logistic map. 

Noise effect on coherent regimes. Let us consider 
populations of noisy and globally coupled identical one- 
dimensional maps. We choose this system because it is 
sufficiently simple to be analytically treated, and, at the 
same time, can be considered as a prototype for inquiring 
new phenomena. 

The fact that the mean field can exhibit complex dy- 
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namics is well known |25|, |2(| |271 |2Jg and in particular 
we will focus on the synchronous regimes appearing for 
sufficiently strong coupling. 

The equation of each population element is defined as 
follows: 



(1-fc) /fa) + *</(*) >+&(*) j = l,2,. 



where Xj € 



.AT, 
(1) 

is the state of the j-th population element, 
/ : M — > R is a smooth function defining the uncoupled 
elements dynamics and (f(x)) := 53jLi f( x j)/N is the 
average over the population. Every map is subject to a 
noise £j(t), that is determined according to a distribution 
of assigned moments (in particular, it does not need to be 
Gaussian). The noise terms are independent and delta 
correlated, that is: 



cascade, crossing several periodic windows and undergo- 
ing a period-halving scenario. This simplification of the 
mean field dynamics does not however reflect on the in- 
dividual elements of the population, whose dynamics is 
more and more smeared out with the increase of the mi- 
croscopic noise. 

As we will later phrase in more rigorous terms, the fact 
that the average has a regular behavior is a consequence 
of the fact that, in the population, a large number of si- 
multaneous realizations of the noise occurs, so that only 
statistical averages of the microscopic stochastic process 
are relevant. 

As a second example of bifurcations induced by micro- 
scopic noise we consider an "excitable map" of the form 
(for other definitions of "excitable map" see Refs. HEIl): 



where a 2 is the variance of the noise terms distribution, 
measuring the intensity of the microscopic noise. 
Keeping all the parameter fixed, we will study the asymp- 
totic behaviour of the mean field X — ( x ) of large pop- 
ulations when cr, that will be from now on our control 
parameter, is changed. 

Figure H shows the bifurcation diagram of X in the case 
of half a million logistic maps in the chaotic regime, to- 
gether with the phase portrait of an individual element 
of the population. When the noise intensity is small, all 
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FIG. 1: Phase portrait of the mean field X (big dots) of 
TV = 2 19 logistic maps of the form f(x) = 1 — ax 2 (a = 1.57, 
k — 0.7) (and of one population element x (small dots)) as a 
function of the noise standard deviation a. The noise terms 
are generated according to a uniform distribution. 

the elements of the population evolve coherently on a 
chaotic attractor. In the limit of zero noise, indeed, all 
the elements are perfectly synchronized and have a com- 
mon chaotic trajectory (the assumption of strong cou- 
pling prevents the formation of clusters) . For larger noise 
intensities, the mean field displays an inverse bifurcation 



f(x) — (a x + 7 x 3 ) e 



-fix 1 



(2) 



The parameters are chosen in a region where the origin 
is the only fixed point, but, if the system is initiated 
far enough from it, there is a chaotic transient. The 
effect of noise on the single map is thus of exciting it 
over threshold, so that a complex dynamics takes place. 
Correspondingly, the mean field has, for small noise in- 




FIG. 2: Phase portrait of the mean field X of TV = 524288 
excitable maps of Eq. (a = .4, /3 = 1, 7 = 8, k = 0.9) as a 
function of the noise standard deviation a. The noise terms 
are generated according to a Gaussian distribution. 



tensity, fluctuations above zero scaling as a/y/N, up to a 
critical point, where the mean field starts displaying large 
amplitude chaotic oscillations (Figure |5J • The asymp- 
totic dynamics simplifies for higher values of cr, the mean 
field going through a backwards bifurcation cascade up 
to a steady state. For very large noise values, the fixed 
point drops again to zero, as a consequence of the fact 
that the map in Eq. J5J) is odd and a strong noise causes 
the population to spread symmetrically around the mean 
field. 
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Order parameter expansion. Let us now address the 
problem from a mathematical viewpoint, performing a 
change of variables that allows us to decouple the macro- 
scopic effect of noise from the dynamic "skeleton" fur- 
nished by the uncoupled map equation. This is achieved 
expressing the position of each population element in 
terms of the mean field and of its displacement from it: 



xj = X 



J = 1,2,. 



,N. 



(3) 



We can now substitute Eq. into the uncoupled ele- 
ment equation and expand it in series around the mean 
field, thus obtaining: 



f( x j 



OO ^ 

/po+E i^/ws' 



(4) 



9=1 



where V q f indicates the <?-th derivative of the function 
/. The equation for the mean field can now be obtained 
directly from its definition, thus getting: 



OO 

E 

9=1 



1 



V«f(X) 



(5) 



We have at this point discarded the term ( £ ) that van- 
ishes in the limit of infinite population size and, other- 
wise, plays the role of a macroscopic noise acting on the 
mean field and scaling as 1/ VN, according to the law of 
large numbers. 

Let us now define a set of new order parameters: 



fin 



(6) 



and compute their evolution by making use of: 
1 



(1-fc) 



OO 

E 

P =i 



and of the fact that, being the positions and the noise 
uncorrelated variables, in the limit N — > oo holds: 
( h(X, e) ) = ( h(X, e) ) ( ) . From simple algebra thus 
follows: 



0.„ 



9 



(l-fc)*"*g- 



oo 1 



Lp=i 



(7) 



where m q = ( £ q ) is the g-th moment of the noise distri- 
bution. 

As a first approximation, valid for high coupling 
strength, the mean field dynamics can be described by 
a scalar equation: 



X»f(X) + 



E 

9=1 



1 



V*f(X) m q 



(8) 
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FIG. 3: Phase portrait, as a function of the (uniformly dis- 
tributed) noise standard deviation a, of a) the mean field 
of a population (N = 524288) of logistic maps of the form 
f(x) = 1 — ax 2 (a = 1.57) for high coupling (k = 0.9); b) the 
mean field of the same population for weaker coupling (k = 0.7); 
c) the zeroth-order expansion Eq. ||HJ; the second order expan- 
sion Eq. 1101 . For lower coupling values, the qualitative agree- 
ment improves if further order parameters are taken into ac- 
count. 



which is obtained by Eqs. (J5J and Q in the limit k — ► 1. 
The two left-hand pictures of Figure show the mean 
field dynamics for the population of logistic maps under 
high coupling (k = 0.9, top) and for the effective mean 
field dynamics Eq. (jSJ (bottom), that in this case takes 
the simple form: 



X 



1 



a a' 



aX J 



(9) 



This equation provides a rescaling of the average popula- 
tion dynamics to that of a single logistic map of the form 
b + aX 2 . 

A more difficult case to analyse is that of the excitable 
maps, since in this case the sum in Eq. (jSJ contains an 
infinite number of terms. Any finite truncation thus in- 
troduces a further approximation, whose validity holds 
for sufficiently weak noise. Figure 01 reports the mean 
field bifurcation diagram for the population (left) and for 
the sixth-order truncation of Eq. JSJ (right), that is able 
to describe the onset of the collective chaotic oscillations 
up to the first periodic windows, ft is worth noticing 
that the terms of high order in the series of Eq. © are 
small for low noise intensity, so being the noise distribu- 
tion moments. Therefore, increasing the noise variance 
from zero causes the nonlincarities of the uncoupled map 
to become progressively important in the mean field dy- 
namics description. In general, we can sec that Eq. (jSJ) 
accounts for the influence of the noise distribution fea- 
tures (described by its moments) on the dynamics of an 
uncoupled element. If this is a polynomial, the g-th mo- 
ment will affect the coefficients of the terms of order less 
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FIG. 4: Phase portrait as a function of the (normally dis- 
tributed) noise standard deviation a of the mean field of a popu- 
lation (N = 2*9) of excitable maps described by Eq. (a = 6, 
(3 = 1, 7 = 8) for k = 0.9 (top left) and for Eq. JSJ truncated at 
the sixth order. Lower order truncations lead to less accurate 
approximations (not shown). 

than q. It is however important to stress that no require- 
ments have been made on the specific shape of the noise 
distribution, so that distributions having the same low- 
order moments have similar macroscopic effects for weak 
noise. 

Let us now consider the case of lower coupling and 
include in the effective description also the terms in (1 — 
k) 2 (those in (1 — k) cancel out for symmetrical noise 
distributions, like those considered so far): 

x^ f{x) + V - vf(x) n q 

tt q h- m q + (i - kf m q _ 2 (10) 

P =i 1 ■ 

If / is a polynomial of order n, the reduced system Eq. 
(|10fl is a closed system of 2n + 1 variables, which thus 
determine the dimensionality of the mean field dynamics. 
Again, when / is not a polynomial, any truncation of 
Eq. (|10|) introduces errors which become bigger when the 
noise intensity increases. Nevertheless, the moments of 
high order become negligible when a is small and the 
mean field dynamics can still be described in the weak 
noise regimes. 

Going back to the previously considered population of 
logistic maps, we lower the coupling strength. Equation 
(|10fl takes now the following form: 

Xh l-aX 2 -an 2 

fl 2 ^ a 2 + {l-k) 2 a 2 (AX 2 Vl 2 -ti 2 2 + Vl i ) (11) 
Q 4 ^ m 4 + 6(1 -k) 2 a 2 a 2 (AX 2 VL 2 - f2|+0 4 ) 

The two right-hand images of Figure|3]compare the mean 
field dynamics for k = 0.7 (top) with that of the second 
order approximation Eq. i|ll|) (bottom), the mean field 
dynamics can be directly compared with the zeroth-order 
approximation (bottom left of the same figure), that is in- 
dependent of k. It is evident that the accordance with the 



population bifurcation scenario is improved when more 
order parameters are introduced in the description. In 
particular, Eq. l(TT|l reproduces the shift of the period- 
doubling cascade toward lower values of the noise inten- 
sity, taking place when the coupling is weakened. 

Conclusions. In this Letter we have discussed the phe- 
nomenon of macroscopic bifurcations induced by micro- 
scopic additive noise in large populations of globally and 
strongly coupled maps. We have shown that the macro- 
scopic dynamics changes with the noise intensity and that 
in the case of high coupling the macroscopic bifurcations 
can be reproduced by a low-dimensional map, where the 
mean field is coupled to some additional order param- 
eters. This order parameter reduction accounts for the 
effects of noise on the dynamical "skeleton" given by the 
uncoupled element equation. The proposed method holds 
for any smooth map and is largely independent from the 
specific characteristics of the noise distribution, allowing 
to understand how the noise moments interact with the 
nonlinearities of the uncoupled map equation. As exam- 
ples of application, we have considered logistic maps in 
the chaotic regime and excitable maps, and we have ad- 
dressed the validity of different approximation schemes. 
Although populations of scalar maps have been consid- 
ered in order to simplify both the analytical and the 
numerical work, studies on continuous-time systems (in 
particular, Refs. [ill together with our preliminary 
results indicate that microscopic noise might have a sim- 
ilar effect on strongly coupled continuous time systems, 
inducing macroscopic bifurcations of the mean field. The 
method presented here could be combined with the order 
parameter expansion proposed in Ref. |30| for investigat- 
ing the different roles of microscopic noise and intrinsic 
parameter diversity on populations of dynamical systems. 
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